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We study the properties of the light pseudoscalar mesons in a three flavor chiral quark model 
with nonlocal separable interactions. We concentrate on the evaluation of meson masses and decay 
constants, considering both the cases of Gaussian and Lorentzian nonlocal regulators. The results 
are found to be in quite good agreement with the empirical values, in particular in the case of 
the ratio fx/ /-tt and the anomalous decay ty° 77. In addition, the model leads to a reasonable 
description of the observed phenomenology in the 77 — 77' sector, even though it implies the existence 
■^|- ■ of two significantly different state mixing angles. 
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The properties of the light pseudoscalar mesons (i.e. the pions, kaons and etas) provide a suitable ground for the 
study of the basic nonperturbative features of Quantum Chromodynamics (QCD). As well known, the QCD Lagrangian 
shows an approximate U(3)l (£> U(3)r chiral symmetry, which is spontaneously broken down to U(3)y in the low 
momentum, nonperturbative regime. The fact that, instead of nine, only eight pseudoscalar quasi-Goldstone bosons 
are observed in Nature is usually explained in terms of the so-called U(\)a anomaly. This anomaly is again related to 
nonperturbative aspects of QCD, and it is believed to be mainly responsible for the rather large 77' mass. Unfortunately, 
so far it has not been possible to obtain detailed information about the properties of the light pseudoscalar mesons 
" l ■ directly from QCD, and most of the present theoretical work on the subject relies in low energy effective theories. 
' Among them the Nambu-Jona-Lasinio (NJL) model [1] and its three flavor extensions [2, 3, 4, 5] are some of the 
most popular ones. In the NJL model the quark fields interact via local effective vertices which are subject to 
chiral symmetry. If such interaction is strong enough chiral symmetry is spontaneously broken and pseudoscalar 
Goldstone bosons appear [6]. As an improvement of the local NJL scheme, some covariant nonlocal extensions 
have been studied in the last few years [7]. Nonlocality arises naturally in the context of several quite successful 
approaches to low-energy quark dynamics as, for example, the instanton liquid model [8] and the Schwinger-Dyson 
resummation techniques [9] . It has been also argued that nonlocal covariant extensions of the NJL model have several 
advantages over the local scheme. Indeed, nonlocal interactions regularize the model in such a way that anomalies are 
• . preserved [10] and charges properly quantized, the effective interaction is finite to all orders in the loop expansion and 
there is not need to introduce extra cut-offs, soft regulators such as Gaussian functions lead to small next-to-leading 
order corrections [11], etc. In addition, it has been shown [12, 13] that a proper choice of the nonlocal regulator and 
the model parameters can lead to some form of quark confinement, in the sense that the effective quark propagator 
has no poles at real energies. 

Until now, most of the research work on nonlocal chiral models has been restricted to the flavor SU(2) sector 
including applications to the baryonic sector [14, 15] and to the study of the phase transitions at finite temperature 
and densities [16]. The aim of the present paper is to extend this type of models so as to include strange degrees of 
freedom, and to analyze the predictions for the masses and decay constants for the pions, kaons and the 77 — 77' system. 

This article is organized as follows. In Sec. II we present the general formalism and derive the expressions needed 
to evaluate the different meson properties. The numerical results for some specific nonlocal regulators together with 
the corresponding discussions are given in Sec. Ill, while in Sec. IV we present our conclusions. Finally, we include 
an Appendix with some details concerning the evaluation of quark loop integrals. 
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II. THE FORMALISM 
A. Effective action 



We start by the Euclidean quark effective action 



Se = J d A x jv>(x) [-i-f^ + m c ] ip{x) - ^ [jf (x) jf (x) + jf (x) jf (x) 



T A afcc [jf (x)jf (x)jf(x) - 3 jf(x)jf(x)j c p (x) 



(1) 



where V' is a chiral £/(3) vector that includes the light quark fields, ip = (u d s) T , while rh c — diag(m u , m<j, m s ) is 
the current quark mass matrix. We will work from now on in the isospin symmetry limit, in which m u — nid- The 
currents jf' p (x) are given by 



jf (x) = J d*y d*z r(y - x) r(x - z) tp(y) X a i>(z) , 
jf (x) = y d 4 y d 4 z r(y - x) r(x - z) $(y) i-f 5 X a ip(z) , 



where the regulator r(x — y) is local in momentum space, namely 

d 4 p 



r( 



x-y) = J 



i(x-y)p 



r(p) 



(2) 
(3) 

(4) 



and the matrices A a , with a = 0, ..,8, are the usual eight Gell-Mann 3x3 matrices — generators of SU(3) — plus 
A = \/2/3 11-3x3 • Finally, the constants A abc are defined by 



A a bc — g| Ejjfc £mn! (A a ) im (Ati)j n (A c ) fe; . 



(5) 



The corresponding partition function Z = j DipDip exp[— Se] can be bosonized in the usual way introducing the 
scalar and pseudoscalar meson fields a a (x) and 7r a (x) respectively, together with auxiliary fields S a (x) and P a {x). 
Integrating out the quark fields we get 



Z = j Va a Vn a dot A J VS a VP a cxp j d 4 x (a a S a + TT a Pa) 



exp 



G H 

— (S a S a + P a Pa) + —^A a b c (S a SbS c — 3S a PbPc 



}' 



where the operator A reads, in momentum space, 

A(p,p') = {-j> + m c )(2iT)U^{p-p l ) + r{p) [a a (p - p') + z 75 ir a (P - p')} X a r(p') . 



(6) 



(7) 



To perform the integration over the fields S a and P a we use the Stationary Phase Approximation (SPA). This 
means to replace the integral over S a and P a by the integrand evaluated at its minimizing values 5 a (<7f,(a;), 7r c (x)) and 
fa(o&(x), 7r c (x)). The latter are required to satisfy 



3H r 

a a + GS a + — A abc S b S c - Pb P c 

n a + GP a -^- A abc S b P c = Q. 



= 



Thus, within the SPA the bosonized effective action reads 

S E = - In dot A - j d 4 x a a S a + n a P a + ^ (s a S a + P P a ) + ^A abc (s a S b S c - 3 S a P b P c ) 



(8) 



(9) 



3 



At this stage we assume that the a a fields can have nontrivial translational invariant mean field values a a while the 
pseudoscalar field cannot. Thus, we write 



a a {x) = a a + 5a a (x) 
ir a (x) = 5ir a (x) . 



(10) 



Note that due to charge conservation only (7 a= o,3,8 can be different from zero. Moreover, 03 also vanishes in the 
isospin limit. After replacing Eqs. (10) in the bosonized effective action (9) and expanding up to second order in the 
fluctuations 5a a (x) and 5ir a {x) we get 



S E = S™ FA + S q E uad + 



Here the mean field action reads 



qMFA 

y(4) - lNc 



i]Vp Trln[p 2 l 3 x 3 + S 2 b)] -\ 



(27T) 



G 



^ ( &iSi + —SiSi ) + —SuSdS. 



(11) 



(12) 



where for convenience we have changed to a new basis in which Si, with i = u,d,s (or equivalently i — 1,2,3) are 
given by 



1 



1 



S u = \ ~ So + S3 H — ■= S$ , Sd = \ - So — S3 H — ■= Ss , S s = \ - So -\ — 7= 5*8 , 



V3 



and similar definitions hold for Oi in terms of a®, 03 and cr§. In Eq. (12) we have also defined £(p) = 
diag(£„(p), Sd(p), S s (p)), with T<i(p) = rrn + <fi r 2 (p), whereas the mean field values Si are given by Si = Si(<fj,0). 
Note that in the isospin limit a u = ad, thus we have £„(p) = Sd(p). 

In order to deal with the mesonic degrees of freedom, we also introduce a more convenient basis defined by 



1 



4>ij = ( A a0a) y - , 

where <p = a, tt, and the indices i,j run from 1 to 3. For the pseudoscalar fields one has in this way 



(13) 



5-Ki 



( tt 77s n 
y/2 V6 \/3 



K- 



K+ 



\ 



.1^ + 1* +1*L K ° 
V2 V6 V3 

V6 V3 I ij 



(14) 



The second piece of the effective action in Eq. (11) — quadratic in the meson fluctuations — can be written now as 



nquad 



where we have defined 



with 



ij,kl 



1 f d A p r 

2 J (2^)4 [ G tiM^ 6a ij(P) foki{-p) + GT jkl (p) Sm^p) 6tt m (- p ) 
d 4 q r 2 (q+y(q-) [(q+ ■ q~) T (<T)] 



(2tt)4 [(q+^ + ^(q+)}[(q-^ + ^(q-)] 



q ± =q±p/2 



H 



G Su Sjk i — Eifcn tjln S n . 



(15) 

(16) 

(17) 
(18) 
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B. Mean field approximation and chiral condensates 

The mean field values a u and a s can be found by minimizing the action Sg FA . Taking into account Eqs. (8), a 
straightforward exercise leads to the following set of coupled "gap equations" : 

H - 

"u + G S u H — S U S S = 



2 

a s +GS s + ^S 2 = 



H - . ("J 
2 



where 



_ f d±p S f (p) r 2 (p) 

The chiral condensates are given by the vacuum expectation values (uu) = (dd) and ( ss ) . They can be easily 
obtained by performing the variation of Z MFA = exp[— Sg FA ] with respect to the corresponding current quark 
masses. For q — u,s one obtains 



C. Meson masses and quark- meson coupling constants 



From the quadratic effective action S q ^ ad it is possible to derive the scalar and pseudoscalar meson masses as well 
as the quark- meson couplings. In what follows we will consider explicitly only the case of pseudoscalar mesons. The 
corresponding expressions for the scalar sector are completely equivalent, just replacing the upper indices "— " by 
"+" . In terms of physical fields, the contribution of the pseudoscalar mesons to S q ^ ad can be written as 

d A p 



Qquad 



(27T) 



GAP 2 ) [7T°(p) 7T°(-p) + 2n + (p) TT-(-p)] 

+G K (p 2 ) [2K°(p) K°(-p) + 2K+(p) K-(-p] 
+G v {p 2 ) r,(p) V(-P) + G v ,(p 2 ) rj'ip) V '(-p)\ . 



Here, the fields rj and i]' arc related to the U(3)v states r]o and rjs according to 

i] = cos V i] 8 — sin 6 V r]o 
rj = sin 6 V > i] S + cos 9 V > rj Q , 



(22) 



(23) 
(24) 



where the mixing angles 9 V:V > are defined in such a way that there is no 77 — 77' mixing at the level of the quadratic 
action. The functions Gp(p 2 ) introduced in Eq. (22) are given by 



GAP 2 ) 
G K (P 2 ) 

gap 2 ) 

G v ,{p 2 ) 



(g + ^SsA' + c-Ap 2 ) 

(G+§S u )-'+C- s (p 2 ) 



G m (p 2 ) + G m {P 2 ) 



Gss(P 2 ) + Go (P 2 ) 



-\[Gos(P 2 f + 



GssiP 2 ) - G AP 2 ) 



+ \[Gos(P 2 f + 



G m {P 2 ) G m (P 2 ) 



(25) 
(26) 

(27) 

(28) 
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where 



2^ V2 



G 0S (P 2 ) = 



Goo(p 2 ) 



6G - AHS U - 2HS S 
2G 2 - GHS S - H 2 S 2 
H(S S — S u ) 
2G 2 -GHS S - H 2 Sl 
6G + ms u - HS 8 



+ C- u (p 2 ) + 2C-(p 2 ) 



2G 2 -GHS s -H 2 Sl 
The meson masses are obtained by solving the equations 

Gp(-m 2 P ) = 



+ c- u ( P 2 )-c-(p 2 ) 



2C- U {p 2 )+C-{P 2 ) 



(29) 



(30) 



with P — 7r, K, r\ and rj , while the tj and i]' mixing angles, which are in general different from each other, are given 

by 



™ G m (p 2 )-G m (p 2 ) 



p 2 =—m 2 , 



(31) 



Now the meson fields have to be renormalized, so that the residues of the corresponding propagators at the meson 
poles are set equal to one. This means that one should define renormalized fields <p(p) = x ^ 2 <p(p) such that, close 
to the poles, the quadratic effective lagrangian reads 



1 



= o (P 2 + ml) <t>{p) 0(-p) . 



2 

In this way, the wave function renormalization constants Zp are given by 

dG P {p 2 ) 



Z P = 



dp 2 



(32) 



(33) 



with P — 7r, K, rj and 77'. Finally, the meson-quark coupling constants Gp q are given by the original residues of the 
meson propagators at the corresponding poles, 



G 



Pq 



A/2 



(34) 



D. Weak decay constants of pseudoscalar mesons 

By definition, the pseudoscalar meson weak decay constants are given by the matrix elements of the axial currents 
A^(x) between the vacuum and the renormalized one- meson states at the meson pole: 

(0\Al(0)\Mp))=ifabPn. (35) 

For a, b = 1 ... 7, the constants f a b can be written as 8 a b f<j>, with (f> = ir for a = 1, 2, 3 and </> = K for a = 4 to 7. In 
contrast, as occurs with the mass matrix, the decay constants become mixed in the a = 0,8 sector. 

In order to obtain the expression for the axial current, one has to gauge the effective action Se by introducing a set 
of axial gauge fields A^. For a local theory, this gauging procedure can be done simply by performing the replacement 



(36) 



However, since in the present case we are dealing with nonlocal fields, an extra replacement has to be performed in 
the regulator r(x — y) [12, 13, 17]. One has 



r(x — y) — > Pexp 



J X 



75 A a A a Js) 



r(x - y) 



(37) 
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where s represents an arbitrary path that connects x with y. In the present work we will use the so-called "straight 
line path" , which means 



+ a(y ll - x^) , 



(38) 



with < a < 1. Once the gauged effective action is obtained, it is straightforward to get the axial current as the 
derivative of such action with respect to A® (x) evaluated at A^(x) = 0. Then, performing the derivative of the 
resulting expressions with respect to the renormalized meson fields, one can finally identify the corresponding meson 
weak decay constants. After a rather lengthy calculation, we find that the pion and kaon decay constants are given 

by 



U - 4/„„(-m 2 ) Zl'\ 
Ik = 2 [fus(-m 2 K ) + fsu(-m 2 K )] Z 



where 



fijip 2 ) - 2N C 



d 4 q ip-qi) 



1/2 
K ! 



(39) 
(40) 



(2tt)' 



P 



[(qir + ^(qi)] [(<fe) 2 +S 2 (<fa)] 



a = l/2 



+ 4N C 



d*q (p-q) dr(q) 
i 2 dq 2 



(2tt)< 



da 



with 



x \r(q a ) [£i(«+)-mi 



la 



[(g+- g -) + S t (g+)S 3 ( g -)] 
(qa) 2 + mqa) 

-■ q+ (1 -a) p , 
q — a p . 



r(qt) £,(</+) 



(41) 



(42) 



In the case of the r\ — r( system, two decays constants can be defined for each component a = 0,8 of the axial 
current [18]. They can be written in terms of the f a b decay constants and the previously defined mixing angles Q n ^ 
as 



= [fasi-m 2 ) cos 9 V - f a0 (-m 2 ) sin 9 n ] Z]l 2 
f% = [fasi-m 2 ,) sind n , + f a0 (-m 2 ,) cos 9 V ,] z]j, 2 . 

Within our model, the decay constants f a b for a, b — 0, 8 are related to the defined in Eq. (41) by 

fss( P 2 ) = \ [2fss{p 2 )+!uu{p 2 )} 
/00(P 2 ) = \ [yuu{p 2 )+fss{p 2 )} 

hs(p 2 ) = /80(P 2 ) = ^ [fuu(p 2 ) fss(p 2 )} ■ 



(43) 

(44) 

(45) 
(46) 

(47) 



It is clear that both the nondiagonal decay constants /os, /so as well as the mixing angles 9 V and 9 V ' vanish in the 
SU(3) symmetry limit. 



E. Anomalous P — > 77 decays 

To go further with the analysis of light pseudoscalar meson decays, let us evaluate the anomalous decays of n°, r/ 
and 7/ into two photons. In general, the corresponding amplitudes can be written as 

A(P -> 77) = e 2 5P77 e^ a p s*/ e* 2 v k<? k% , (48) 

where P = tt°, 77, 77', and fc i; £j stand for the momenta and polarizations of the outgoing photons respectively. 

In the nonlocal model under consideration the coefficients gp 77 are given by quark loop integrals. Besides the usual 
"triangle" diagram, given by a closed quark loop with one meson and two photon vertices, in the present nonlocal 
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scheme one has a second diagram [13] in which one of the quark-photon vertices arises from the gauge contribution 
to the regulator, see Eq. (37). The sums of both diagrams for 7r°, t] and rj' decays yield 



57r77 = I u {ml) Z X J 2 
1 



9vil — 



3\/3 
1 



[5 J„(m 2 ) -2I s (m 2 v )] cos0„ - V2 [5/„(m^) + J a (m^)] sin0„ 



^y 2 



3^ 

where the loop integrals If(m 2 P ) are given by 



[5 J„(m 2 ,) - 2/ s (m 2 / )] sin0„i + V2 [5/ u (m 2 ,) + / s (m 2 ,)] C os0„< 



7 l/2 



r(q - k 2 ) r(q + ki) 



Y [q* + (<?)] [(q + + T? f {q + h)} [(q - fc 2 ) 2 + V)(q - k 2 )] 



(h ■ q) 



(h ■ q) 



(49) 



(50) 



(notice that for on-shcll photons these integrals are only functions of the scalar product [k\ ■ k 2 ) 1 which in Euclidean 
space is equal to —mp/2). In terms of the parameters gp 17 , the corresponding decay widths are simply given by 



r(P -» 77) = J m 3 P a 2 g 2 P ^ , 



(51) 



where a is the fine structure constant. 



F. Low energy theorems 



Chiral models are expected to satisfy some basic low energy theorems. In this subsection we consider some important 
relations such as the Goldberger-Treiman (GT) and Gell-Mann-Oakes-Renner (GOR), showing explicitly that they 
are indeed verified by the present model. 

We start by the GT relation. Taking the chiral limit — > in f uu (— m 2 ) and dG^/dp 2 \ p 2 = _ m 2 appearing in 
Eqs.(39) and (33) we get 

lim fuu{-ml) = \ a ufi Z~J , (52) 

where here, as in the rest of this subsection, the subindex indicates that the corresponding quantity is evaluated in 
the chiral limit (notice that <7 Uj0 = E Ui o(0)). Replacing this expression in Eq. (39) and taking into account Eq. (34) 
we get 

/tt,o Gnqfl = <Ju,o , (53) 

which is equivalent to the GT relation at the quark level in our model. 

Let us turn to the GOR relation. Expanding C~ u (p 2 ) in Eq. (17) to leading order in m u and p 2 we get 

C^ 2 ) ^-,„ + P 2 Z-l . (54) 

a u,0 cr ufi 

To obtain this result we have used the gap equations (19) and the expression for the chiral condensate given in 
Eq. (21). Now using Eq. (54) together with the equation for the pion mass, 

GA-ml) = , (55) 

and taking into account the GT relation (53), one gets 

m u ( uu + dd ) = - fl fi ml , (56) 



which is the form taken by the well-known GOR relation in the isospin limit. 
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Next we discuss the validity of the Feynman-Hellman theorem for the case of the so-called pion sigma term. This 
theorem implies the relation [19] 



dm 2 
dm u 



( 7T | UU + dd | 7T ) 



(57) 



where covariant normalization, (p'\p) = 2E p (2n) 3 S 3 (p' —p), has been used for the pion field. An expression for the 
left hand side of Eq. (57) can be easily obtained by deriving Eq. (55) with respect to the u quark mass. In fact, 



dG^-ml) _ 8G n (p 2 ) 



dm v 



+ 



dG w (p 2 ) 



dp 2 





f dm l\ 


p 2 =—m2 


V dniu) 



(58) 



thus 



dm\ 
dm,. 



dGAp 2 ) 
dm u 

dG^(p 2 ) 
dp 2 



(59) 



On the other hand, within the path integral formalism, one has 

s 3 s E \j] 



( 7T Q I UU + dd I Tr" ) = 



6j(0) 5n a (p) 5ir b {-p) 



X 



n = j = 
p 2 = -ml 



(60) 



where 



(61) 



S 3 E \j] = S E + Jd 4 x[ u(x)u(x) + d{x)d{x) } j(x) , 

Se being the effective action of the model, given by Eq. (1). From the explicit form of Se it is easy to see that 

S E \j] = S E (m u -» m u + j(x)) , (62) 
therefore, using the bosonized form of the effective action in Eq. (11), with S q E uad given by Eq. (22), we get 

dG^p 2 ) 

1 dG^jp 2 ) 

j = 



( 7T | uu + dd | 7r) = 



dGi\j] 



dj(o) 



P 2 = -ml 



Z n 1 dm v 



dm u 



dG^p 2 ) 
dp 2 



(63) 



p^ — — 



Comparing Eq. (63) with Eq. (59) we see that the FH theorem, as it should, holds in the present model. Moreover, 
using the GOR relation Eq. (56), we obtain up to leading order in m u 



( 7T | UU + dd | 7T ) 



(uu + dd) _ m 2 
flo m u 



(64) 



To conclude, let us analyze in the chiral limit the coupling g^^. Expanding the integrand of I u (m 2 ) in powers of 
ki, fc 2 and taking the limit m 2 — ► one obtains [20] 



7 l/2 
_ ^7r,0 

Now, taking into account the GT relation, one finally has 

1 

which is the expected result according to low energy theorems and Chiral Perturbation Theory. 



(65) 



(66) 
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III. NUMERICAL RESULTS 

In this section we discuss the numerical results obtained within the above described nonlocal model. Our results 
include the values of meson masses, decay constants and mixing angles, as well as the corresponding quark constituent 
masses, quark condensates and quark-meson couplings. The numerical calculations have been carried out for two 
different regulators often used in the literature: the Gaussian regulator 

r(p 2 ) = cxp (-p 2 /2A 2 ) (67) 

and a Lorentzian regulator 

r(p 2 ) = (l+p 2 /A 2 y 1 , (68) 

where A is a free parameter of the model, playing the role of an ultraviolet cut-off momentum scale. Let us recall 
that these regulators are defined in Euclidean momentum space. 

A. Fits to physical observables 

The nonlocal model under consideration includes five free parameters. These are the current quark masses fh and 
m s (rh — m u = m<j), the coupling constants G and H and the cut-off scale A. In our numerical calculations we have 
chosen to fix the value of fh, whereas the remaining four parameters have been determined by requiring that the 
model reproduces correctly the measured values of four physical quantities. The observables we have used here are 
the pion and kaon masses, the pion decay constant f v , and a fourth quantity, chosen to be alternatively the rf mass 
or the squared rj' — ► 77 decay constant, g 2 ,^- In the case of the Gaussian (Lorentzian) regulator, we find that for 
fh above a critical value m crit ~ 8.3 MeV (3.9 MeV) the quark propagators have only complex poles in Minkowski 
space. This can be understood as a sort of quark confinement [12, 13]. In contrast, for m < m cr i t one finds that u and 
d quark Euclidean propagators do have at least two doublets of purely imaginary poles (i.e. real poles in Minkowski 
space). 

In Table I we quote our numerical results for several quantities of interest. Besides the obtained values of meson 
masses and decay constants, we include in this table the corresponding results for quark condensates, constituent quark 
masses and quark-meson coupling constants. We have taken into account both Gaussian and Lorentzian regulators, 
considering in each case values for the light quark mass fh above and below m cr it- Sets GI, GIV, LI and LIV have 
been determined by fitting the free parameters so as to reproduce the empirical value of m,/, while for sets Gil, GUI, 
LII and LIII the 77' mass is obtained as a theoretical prediction and the parameters have been determined by adjusting 
5 2 / 77 to its present central experimental value. The last column of Table I shows the empirical values of masses and 
decay constants to be compared with our predictions. 

As stated in Sect. II. C, the meson masses are obtained by solving the equations G p(—m 2 p ) = for P = it, K , 77 and 
77'. Now, to perform the corresponding numerical calculations, one has to deal with the functions Cij(p 2 ) evaluated 
at Euclidean momentum p 2 — —m 2 P . These functions, defined by Eq. (17), include quark loop integrals that need to 
be treated with special care when the meson mass exceeds a given value mp > 2 Si, where Si is the imaginary part 
(in Euclidean space) of the first pole of the quark propagator. In practice this may happen in the case of the 77' state, 
and physically it corresponds to a situation in which the meson mass is beyond a pseudothrcshold of decay into a 
quark-antiquark state. A detailed discussion on this subject is given in the Appendix. In particular, it is seen that for 
fh > m cr it (i.e. if the propagator has no purely imaginary poles) the quark loop integrals can be regularized in such 
a way that their imaginary parts vanish identically, and consequently the rj' width corresponding to this unphysical 
decay is zero. In contrast, for rh < m cr i t the width is in general nonzero, and the presence of an imaginary part 
implies that the condition G v >(—m 2 ,) = cannot be satisfied. In any case, it is still possible to define the 77' mass 
by looking at the minimum of |G r? /(— p$)\. The situation is illustrated in Fig. 1, where we plot the absolute values 
of the functions Gp(-Pq) for Sets GI and GIV (upper and lower panel, respectively). For Set GI (fh > m cri t) the 
uu pseudothrcshold is reached at about 1 GeV, above the 77' mass, thus the quark loop integrals are well defined and 
G v >(—m 2 ,) = 0. It can be seen that this is also the situation for the parameter Sets Gil, GUI, LII and LIII. On the 
other hand, for Set GIV (fh < m crit , lower panel in Fig. 1) the uu pseudothreshold is reached at ~ 750 MeV, well 
below the r) 1 mass. As stated, in this case one can define m v ' by looking at the minimum of the function 
which is represented by the dashed-dotted curve. Though different from zero, the (unphysical) 77' width is relatively 
small, and the minimum (which, as required, lies at po = 958 MeV) is close to the horizontal axis in the chosen GeV 2 
scale. Therefore we do not expect our results to be spoiled by confinement effects not included in the model. In any 
case, we believe that for fh < m cri t the values of physical parameters related to 77' decay may not be reliable, and 
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they have not been included in Table I. A similar situation occurs for Set LIV. Finally, in the case of Set LI it turns 
out that while m > m crit (no purely imaginary poles), the uu pseudothrcshold occurs below the 77' mass. In this case, 
in order to evaluate the loop integrals we have followed the regularization procedure described in the Appendix, in 
which the imaginary parts of quark loop integrals vanish, leading to a real r( pole. 

By examining Table I it can be seen that, for the chosen values of m, the results for the quark condensate (uu), 
the ratio (ss)/ (uu) and the constituent quark masses S 9 (0) are similar to the values obtained within most theoretical 
studies [6, 9]. The most remarkable differences between the results corresponding to both regulators are found in the 
values of the current quark masses and the quark condensates. For the Gaussian regulator the parameters fh and 
m s are found to be about 8 MeV and 200 MeV respectively, while for the Lorentzian regulator the corresponding 
values are approximately 4 MeV and 100 MeV. Note, however, that the obtained values of the ratio m s /fh are very 
similar and somewhat above the present phenomenological range (m s /m u ) emp = 17 — 22 [18]. On the contrary, as 
expected from the GT relation and its generalizations, quark condensates are found to be higher for the Lorentzian 
regulator parameter sets. In any case, the results for both regulators are in reasonable agreement with standard 
phenomenological values [21] and the most recent lattice QCD estimates [22]. 

The predicted values for the kaon decay constant are also phenomenologically acceptable. Indeed, the prediction 
for the ratio fx/ fir turns out to be significantly better than that obtained in the standard NJL, where the kaon and 
pion decay constants are found to be approximately equal to each other [6] in contrast with experimental evidence. 
It is also worth to notice that we obtain a very good prediction for the 7r° — > 77 decay rate. In this sense the nonlocal 
model shows a further degree of consistency in comparison with the standard local NJL model, where the quark 
momenta in the anomalous diagrams should go beyond the cutoff limit in order to get a good agreement with the 
experimental value [7] . 

Our results for 77 and 77' masses and 77, 77' — > 77 anomalous decays require a more detailed discussion. In the Gaussian 
regulator case, it is seen that Set GI, while fitting m v i, leads to a rather large value for g^/ 77 . On the contrary, for set 
Gil, which fits the value of <7^ 77 , the rj' mass decreases up to a value of about 880 MeV (of course, one can also choose 
intermediate sets between this two). In addition, it is found that the fit to <7~/ 77 has a second solution, namely the 
parameter Set GUI, which leads to a 77' mass of about 1 GeV. However, in this case m,< is found to be very close to 
the pseudothreshold point, and as a consequence both the values of the r/-quark coupling G v ' q and the decay constant 
Sy 77 are quite sensible to small changes in the parameters. On the other hand, for all four Gaussian regulator sets 
the results for m n and g nll do not change significantly. The values for the r\ mass are found to be in relatively good 
agreement with experiment, while fl^ 77 turns out to be somewhat large. In the case of the Lorentzian regulator, the 
above described situation becomes strengthened: set LI leads to an unacceptably large value for g v < 71 (notice that, as 
stated, here m,/ is above the pseudothrcshold), while set LII leads to a too low 77' mass. Set LIII seems to reproduce 
reasonably all measurable quantities, but as in the Gaussian regulator case, the result for g,/ 77 is highly unstable. In 
general, as a conclusion, one could say that the Gaussian regulator is preferred, leading to a reasonable global fit of 
all quantities considered here. 

It is worth to mention that the chosen value of m cannot be very far from the values considered in Table I. For 
higher m one would obtain too low values for the quark condensates. On the other hand, if m > m crit the 77' mass 
turns out to be very far above the uu pseudothreshold implying the existence of a large unphysical 77' width. 



The problem of defining and (indirectly) fitting mixing angles and decay constants for the 77 — 77' system has been 
revisited several times in the literature (see Ref. [23], and references therein). On general grounds one has to deal 
with two different state mixing angles 9p and four decay constants f P , where P = 77,77' and a = 0,8. Standard 
analyses, however, used to parameterize the mixing between both meson states and decay constants using a single 
parameter (a mixing angle usually called 9), and introducing two decay constants, fs and /o, related to .g,, 77 and 
<jy 77 through low energy equations analogous to Eq. (66). In the last few years, the analysis has been improved 
(mainly in the framework of next-to-leading order Chiral Perturbation Theory), and several authors have considered 
the possible non-orthogonality of (/®,/®/) and (/",/"/) [24, 25], as well as that of the states 77 and 77' [26]. For the 
sake of comparison, we follow here Ref. [24] and express the four decays constants fp in terms of two decay constants 
f a and two mixing angles a , where a = 0,8. Namely, 



In our model, the decay constants fp can be calculated from Eqs. (43-47). As shown below, it turns out that the 
angles 9s and 9q are in general different. This also happens with the mixing angles 9 V and 9 V >, whose expressions are 



B. 77 and 77' mixing angles and decay constants 




(69) 
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given in Eq. (31). Notice that these are consequences of the (rather strong) p 2 dependence of the functions C\j(p 2 ) 
and fij(p 2 ) defined in Eqs. (17) and (41), respectively. 

Our numerical results for the the parameters f a , 9 a introduced in Eq. (69) and the mixing angles 9 n . v > are collected 
in Table II. In the last column of the Table we quote the ranges in which the parameters f a , 9 a fall within most 
popular theoretical approaches. We have taken these values from the analysis in Ref. [23], in which the results from 
different parameterizations have been translated to the four-parameter decay constant scheme given by Eq. (69). By 
looking at Tabic II it is seen that the results corresponding to Gaussian regulator sets lie within the range of values 
quoted in the literature, while in the case of the Lorentzian regulator the most remarkable difference corresponds to 
Set LI, which leads to a large value of /o (this can be related with the unacceptably large value of g,,' 77 discussed 
above). The mixing angles 9 a can be compared to those obtained in a recent Bethe-Salpeter approach calculation [27], 
which leads to a somewhat larger (absolute) value of 9 . As stated, for Sets GIV and LIV the values for the 77' decay 
constants have not been computed in view of the unphysical finite width acquired by the 77' meson. 

Finally, notice that for most parameter sets (both Gaussian and Lorentzian) the mixing angle 9 V gets a small 
positive value, whereas the angle 9 r y lies around —45°. That is to say, we find our mixing scheme to be very far from 
the approximation of a single mixing angle 8. For this reason we believe there is no reason to expect 9 V to be within 
the usually quoted range 8 = -(10° - 20°). 

IV. CONCLUSIONS 

In this work we have studied the properties of light pseudoscalar mesons in a three flavor chiral quark model 
with nonlocal separable interactions, in which the U(1)a breaking is incorporated through a nonlocal dimension nine 
operator of the type suggested by 't Hooft. We consider the situation in which the Minkowski quark propagator has 
poles at real energies as well as the case where only complex poles appear, which has been proposed in the literature as 
a realization of confinement. We concentrate on the evaluation of the masses and decay constants of the pseudoscalar 
mesons for two different nonlocal regulators, namely Gaussian and Lorentzian. 

As general conclusions, it is found that in this model the prediction for the ratio fx/fir turns out to be significantly 
better than that obtained in the standard NJL, where the kaon and pion decay constants are found to be approximately 
equal to each other in contrast with the experimental evidence. In addition, the model overcomes the standard local 
NJL problem of treating the anomalous quark loop integrals in a consistent way. With respect to the 77 — 77' system, 
according to our analysis the Gaussian regulator seems to be more adequate than the Lorentzian one to reproduce the 
observed phenomenology. In general, our results are in reasonable agreement with experimentally measured values, 
and the global fits can be still improved by considering regulators of more sophisticated shapes. Alternatively one 
might consider adding degrees of freedom not explicitly included in the present calculations such as explicit vector 
and axial-vector interactions, two-gluon components for the 77 and 77' mesons, etc. Finally, it is worth to remark that 
our fits lead to an 77 — 77' system in which the U (3) states 77s and 770 are mixed by two angles 8 n and 9 V > that appear 
to be significantly different from each other. 
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Appendix: Evaluation of quark loop integrals 



In this Appendix we describe some details concerning the evaluation of the quark loop integrals dj(p 2 ) defined in 
Eq. (17). Notice that in the nonlocal chiral model analyzed in this work all four-momenta are defined in Euclidean 
space. However, in order to determine the meson masses, the external momenta p in the loop integrals have to be 
extended to the space-like region. Hence, without loss of generality, we choose p = (0, ipo), and use three-momentum 
rotational invariance to write the quark loop integrals as 

Cij{ ~ Po) = J dq3 dqi 93 + [( q -r + m q -)] ' (A - 1} 

where (73 = \q\. The explicit form of -FVjfe, <74,Po), which depends on the meson under consideration, can be easily 
obtained by comparing Eqs. (A.l) and (17). In principle, the integration in Eq. (A.l) has to be performed over the 
half-plane 93 e [0, 00), 94 e (—00, 00). For sufficiently small values of po the denominator does not vanish at any point 
of this integration region. However, when po increases it might happen that some of the poles of the integrand pinch 
the integration region, making the loop integral divergent. In such cases one has to find a way to redefine the integral 
in order to obtain a finite result. In practice, we need to extend the calculation of the loop integrals Cy(— p^) to 
relatively large values of po only when trying to determine the 7/ mass. For this particle, one only has to deal with 
integrals in which i = j, therefore we will restrict here to this case and the indices i,j will be dropped from now on. 

Let us start by analyzing the zeros of the denominator in the integrand in Eq. (A.l). This denominator can be 
written as D = D + D~ , where 

D ± = {q ± f + T.\q ± ) . (A.2) 

It is clear that the zeros of D are closely related to the poles of the quark propagator S(q) = [— $ + ^(q)]^ 1 . In what 
follows we will assume that the regulator is such that the propagator has a numerable set of poles in the complex q 
plane, and that there are no cuts. It is not hard to see that these poles appear in multiplets that can be characterized 
by two real numbers (S", S"), with > 0, S" > 0. The index z/eN has been introduced to label the multiplets, with 
the convention that v increases for increasing S\ . It is convenient to distinguish between two different situations: 
(a) there are some purely imaginary poles, i.e. there are one or more v for which S" = 0; (b) no purely imaginary 
pole exists, i.e. S" > for all v. It can be shown that purely imaginary poles show up as doublets located at 
Euclidean momentum (yfq 2 )" = ±iSf, while complex poles [28] appear as quartets located at {\fcf) v = S" ± i S" 
and (\fq 2 ) v = —S"±i S" . Clearly, the number and position of the poles depend on the specific shape of the regulator. 
For the Gaussian interaction, three different situations might occur. For values of a below a certain critical value a c , 
two pairs of purely imaginary simple poles and an infinite set of quartets of complex simple poles appear. It is possible 
to check that in this case one of the purely imaginary doublets is the multiplet which has the smallest imaginary part 
(y = 1, according to our convention). At a = a c , the two pairs of purely imaginary simple poles turn into a doublet of 
double poles with S r = 0, while for a > a c only an infinite set of quartets of complex simple poles is obtained. In the 
case of the Lorentzian interactions, there is also a critical value above which purely imaginary poles at low momenta 
cease to exist. However, for this family of regulators the total number of poles is always finite. 

As stated, for low enough external momentum po the integrand in Eq. (A.l) does not diverge along the integration 
region. As p increases, the first set of poles to be met is that with the lowest value of Si, namely (\/<f) u=1 ■ In the 
calculation of the meson properties mentioned in the main text, we deal with relatively low external momenta, so that 
the effect of higher poles is never observed. Thus, in order to simplify the discussion we will only consider in what 
follows the first pole multiplet, dropping the upper index v. The extension to the case in which other sets of poles 
become relevant will be briefly commented at the end of this Appendix. 

The denominator D vanishes when D + = and/or D~ = 0, i.e. when 

(q + f = ql + ql ~ f + TO = S 2 r - Sf ± 2*^5, (A.3) 



and/ or 



2 

(q~ f = ql + ql - ^ - mpo = s 2 r - sf ± 2iS r s t . (A.4) 



Solving these equations for 94 we get in general eight different solutions. Four of them are given by 
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where 

and the other four solutions arc = — Rx^g^ + zlm^'), with z = 1 . . . 4. In Eqs. (A. 5) and (A. 6) g^ 3,1 ^ correspond 

to the zeros of D + and g^ 4,2 ' to those of D~ , while a similar correspondence holds for . For purely imaginary 
poles one has S r = 0, hence only four independent solutions exist. 

If po is relatively small, the distribution of the poles in the complex g4 plane is that represented in Fig. 2. Fig. 
2a holds for situation (a), in which the poles in the first multiplct are purely imaginary while Fig. 2b corresponds 
to case (b), in which these poles are complex. In both figures the dots indicate the zeros of D + and the squares 

those of D~ . As we see, for small values of po half of the poles of D + (namely for case (a) and q^ and q^ 

for case (b)) are placed below the real axis, whereas the other half (g^ for case (a) and and qf^ for case (b)) 

lie above it. Something similar happens for the poles of D~ . Now, as p increases, the poles move in the direction 

(2) (3) 

indicated by the arrows. For a certain value of po, the poles q\ and q\ meet on the real P4 axis (the same obviously 

happens with g^ and g^ in case (b)), thus pinching the integration region of the (qs, 94) integral in Eq. (A.l). The 
location of this so-called "pinch point" in the (g3,g4) plane, which we denote by (q^qf), is given by the solution of 
Irngf = Imgf } = 0: 

(«s>«*)- ^ ^ ,±—j ■ (A.8) 

As we see, for po < 2Si there is no pinch point (actually, it occurs for a complex value of g3, outside the integration 
region in Eq. (A.l)), while for p > 2Si one or two pinch points exist depending on whether S r = (case (a)) or 
S r 7^ (case (b)). In this way, for p > 2 Si the integral in Eq. (A.l) turns out to be ill-defined. 

In order to find a proper regularization procedure, let us analyze a simpler situation in which the problem might 
be solved in Minkowski space through the usual "i e" prescription. We consider the loop integral that appears in the 
usual Nambu-Jona-Lasinio model with three-momentum cut-off, 

I2{ - Pl) =2i ^T ^mr-rn^iel[ { q M r-^ + ie] ' (A ' 9) 

where we have added the subscript M to stress the fact that here the momenta arc defined in Minkowski space. For 
sufficiently small values of po, even in the limit e — > + , the integral is convergent and no regularization is needed. 
Thus one can simply perform the Wick rotation P4 = ipo and take the limit e — > + even before performing the 
integration. One gets in this way 

h ^ = 2 r f?%\< +^ 2TT7 -^2^ 21 ' (A.10) 
J ( 27T r [(<Zb) 2 + m 2 \ [(q E ) + m \ 

which is an integral of the type given in Eq. (A.l). Note that the poles of the propagators are such that this situation 
belongs to case (a), with Si = m. For po > 2 m the straightforward transformation from Minkowski to Euclidean space 
mentioned above cannot be done, since some poles go through the integration contours. The question is whether the 
result of the well-defined Miskowskian integral (A.9) can be still recovered if one starts with the Euclidean integral 
(A. 10), which is ill-defined for p > 2m due to the presence of a pinch point at (^,04) = {yPo/4: — m 2 , 0). It is 
not hard to prove that the answer is yes, once the g4 integration contours and the pole positions are conveniently 
modified. The procedure requires to introduce two small parameters, e and 5, and take the limit 5 — > + , e — > + at 
the end of the calculation. The parameter e is used to shift the poles of D + and D~ (see Fig. 3), whereas 5 is used 
to split the g3 integration interval in three subintervals: the first region corresponds to g3 > g| + 6, the second to 
gf — S < g3 < gf + S and the third to g3 < gf — 6. For each g 3 region we define a different g4 integration contour, as 
represented in Fig. 3 (in the second region, Fig. 3b, also an arbitrary constant k > 1 is introduced). In fact, in the 
first and third regions the limit e — > + can be taken even before performing the integrations. These two g3 regions 
give the full contribution to the real part of the result. However, more care has to be taken with the intermediate 
region, which is responsible for the full contribution to the imaginary part. For example — as it is well known from 
the "i e" Minkowskian formulation — , changing the sign of e does not affect the real part of the result but does change 
the sign of the imaginary part. 
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The prescription just described can be now applied to regularize any loop integral of the form given in Eq. (A.l) 
in which p > 2S"j. Let us consider first the case (a), for which the lowest set of poles has S r = and the pinch point 
is located at ((73,(74) = (\/Po/4 — Sf, 0). Defining C^(-p^) as an integral of the form given by Eq. (A.l) for which 
only one set of purely imaginary poles contributes, we get 



Re 



c^(-pI) 




q%-6 



dq3 



dq 4 



<?i F{q 3 ,q 4 ,p ) 



[((7+) 2 + £ 2 ((7+)] [((7-) 2 + £ 2 (<7-)] 



<?I F(q 3 ,q 4 ,p ) 



Im 



C (a) (-Po) 



^ gg F(gg,0,po) 
2po 



1 + dq 2 



(A.11) 
(A.12) 



Here p$, 8) is the so-called "residue contribution", responsible for the cancellation of the divergence appearing 

in the integrals in (A. 11) in the limit 8 — > + . Its explicit expression reads 



Jo 



d<73 



Re 



93 ^(<?3,<Z4,Po) 



9S 2 (g-) 



[(«+)* + £*(g+)] 



(A.13) 



94=«< 2) ^(^-yif+sF) 



For case (b) we have to extend the previous analysis to the situation in which the poles are complex even in the limit 
e — > + . In this case there is an ambiguity on how to extend the "i e" prescription already in Minkowski space. Here 
we will follow the suggestion made in Ref. [29], in which opposite signs of e are used for each pole and its hermitian 
conjugate (both defined in Minkowski space). In our Euclidean framework, this corresponds to choose different signs 



of e for sets pf 1 



and p)± . It is not hard to see that with this prescription the contributions to the imaginary part 
of the quark loop integral coming from both sets of poles cancel each other. In this way, defining C^ b '(— p^) as an 
integral of the type given in Eq. (A.l) for which only one set of complex poles contributes, we get 



Im 



C^\-pl)]=0 



(A.14) 



For the real part one obtains an expression similar to Eq. (A. 11), just replacing R( a \— p^, 8) by R^ b \— p^, 8), with 



R {b \-p 2 ,8) = 4n r dq 3 Re 
Jo 



<?3 F (l3,H,Po) 



d(q-y 2 



[((7+) 2 +E 2 (^ 



(% + f) 



(A.15) 



(2) 
<74=<?4 



In principle, the extension of the present analysis to the situation in which further sets of poles are relevant is rather 
straightforward. However, some care has to be taken if one has more than one set of purely imaginary poles, since in 
that case double poles might show up for p > Sf + Sf. 

One is faced with a similar problem in the calculation of the rj' — > 77 decay constant g v > 17 , where the loop integral 
Iu(mt)>) defined by Eq. (50) is divergent for m,/ > 2 Si. Though the situation is slightly more involved (one finds 
two pinch points instead of one), our regularization prescriptions can be trivially extended to include this case. In 
the same way, the method could be extended to more complicated situations, as e.g. those found in Schwinger-Dyson 
type of calculations [9, 30]. Finally, we note that although our prescription has some similarities with that used in 
Ref. [13], the regularization procedure is not exactly the same. In the ((73,(74) integral which appears in Eq. (A. 11) 
the excluded region around the pinch point is a slide of infinite size in the (74 direction and size 28 in the (73 direction, 
as opposed to the circular region of radius 8 used in Ref. [13]. This leads to some minor differences in the numerical 
values of the regulated integrals. 
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Gaussian 






Lorentzian 






Set 




GI 


Gil 


GUI 


GIV 


LI 


LII 


LIII 


LIV 


Empirical 


rh 


[ MeV ] 


8.5 


8.5 


8.5 


7.5 


4.0 


4.0 


4.0 


3.5 


(3.4 - 7.4) 


m 3 


[ MeV ] 


223 


223 


223 


199 


112 


110 


112 


100 


(108 - 209) 


A 


[ MeV ] 


709 


709 


709 


768 


1013 


1013 


1013 


1110 




GA 2 




10.99 


11.44 


10.80 


10.43 


14.68 


16.76 


14.96 


14.05 




-HA 5 




295.3 


275.4 


303.7 


305.1 


743.4 


573.8 


720.4 


821.0 




-< uu > 1/3 


[ MeV ] 


211 


211 


211 


220 


275 


275 


275 


288 




-< ss 


[ MeV ] 


186 


187 


185 


204 


297 


307 


299 


314 




S„(0) 


[ MeV ] 


313 


313 


313 


295 


299 


299 


300 


281 




E s (0) 


[ MeV ] 


650 


662 


645 


607 


562 


615 


569 


518 






[ MeV ] 


139* 


139* 


139* 


139* 


139* 


139* 


139* 


139* 


139 


mK 


[ MeV ] 


495* 


495* 


495* 


495* 


495* 


495* 


495* 


495* 


495 


m v 


[ MeV ] 


517 


505 


521 


522 


543 


513 


540 


545 


547 


m v i 


[ MeV ] 


958* 


879 


1007 


958* 


958* 


778 


908 


958* 


958 


Gjrq 




3.28 


3.28 


3.28 


3.09 


3.13 


3.13 


3.13 


2.94 




GKq 




3.47 


3.52 


3.45 


3.21 


3.05 


3.24 


3.08 


2.80 








3.07 


3.03 


3.08 


2.83 


2.74 


2.69 


2.74 


2.49 




G i 




1.62 


2.01 


1.21 




1.36 


2.21 


1.13 








[ MeV ] 


93.3* 


93.3* 


93.3* 


93.3* 


93.3* 


93.3* 


93.3* 


93.3* 


93.3 


Ik/ fir 




1.29 


1.29 


1.29 


1.29 


1.25 


1.28 


1.26 


1.25 


1.22 


2 

9ivyy 


[ GeV~ 2 ] 


0.073 


0.073 


0.073 


0.074 


0.074 


0.074 


0.074 


0.074 


0.075 ± 0.005 


2 

Qiyy-y 


[ GeV- 2 ] 


0.095 


0.106 


0.091 


0.094 


0.072 


0.108 


0.075 


0.075 


0.067 ± 0.006 


2 

9r)' 77 


[ GcV~ 2 ] 


0.141 


0.116* 


0.116* 




0.278 


0.116* 


0.116* 




0.116 ±0.005 



TABLE I: Numerical results for quark effective masses and condensates, and pseudoscalar meson masses and decay parameters. 
The light quark mass rh has been taken as input, and the parameters m s , A, G and H have been chosen so as to reproduce 
the empirical values of pion and kaon masses, the pion decay constant f w and, alternatively, the r\ mass or the measured value 
of g |2 / 77 (marked with *). 



Set 


GI 


Gil 


giii 


GIV 


LI 


LII 


LIII 


LIV 


Theory & 
Phenomenology 


hlh 


1.33 


1.34 


1.32 




1.32 


1.36 


1.31 




(1.2 - 1.4) 


fo/fn 


1.28 


1.18 


1.18 




1.64 


1.18 


1.05 




(1.0 - 1.3) 


0s 


-14.2° 


-19.7° 


-10.0° 




-11.7° 


-25.4° 


-10.6° 




-(22° - 19°) 


0o 


-2.14° 


-5.25° 


-1.17° 




1.66° 


-6.52° 


1.60° 




-(10° - 0°) 


9 V 


4.65 


1.81 


5.69 


5.32 


7.69 


-0.72 


6.90 


8.00 




6 v i 


-50.0 


-47.1 


-52.4 




-48.4 


-44.6 


-50.6 







TABLE II: Numerical results for r) and rf weak decay parameters and state mixing angles, according to the definitions given 
by Eqs. (31) and (69). 
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200 400 600 800 1000 

p [MeV] 



FIG. 1: Absolute values of the inverse meson propagators as functions of the momentum. The full line corresponds to G w , the 
dashed line to Gk, the dotted line to G v and the dashed-dotted line to G v i . The upper panel corresponds to Set GI and the 
lower one to Set GIV. In both cases the insertions show a detail of the low momentum region where the zeros corresponding to 
the ground state pion, kaon and eta mesons occur. 
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FIG. 2: Schematic distribution in the complex t/4-plane of the poles corresponding to the lowest Si set for: a) a < a c ; b) 
a > a c . In both cases the dots indicate the poles of D + and the squares those of D~ . These distributions correspond to a 
value of po < 2Si . The arrows indicate the movement of the poles as po increases. 
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FIG. 3, Integration paths in the complex „ plane for: ,) „ > + S; t» yfi 

qs < \J ?f — m 2 — 8. The constant k is an arbitrary real number satisfying k > 1. 
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